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The "marginal" distribution for measurable coordinate and spin projection is introduced. Then, 
the analog of the Pauli equation for spin-^ particle is obtained for such probability distributions 
instead of the usual wave functions. That allows a classical-like approach to quantum mechanics. 
Some illuminating examples are presented. 
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o 
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I. INTRODUCTION 

- 

Since the early days of quantum mechanics, we have been forced to coexist with complex probability amplitudes 
without worrying about their lack of any reasonable physical meaning. One should not ignore, however, that the 
wave-like properties of quantum objects still raise conceptual problems on whose solutions, a general consensus is far 
from having been reached 

A possible way out of this difficulty has been implicitly suggested by Feynmann who has shown that, by dropping 
the assumption that the probability for an event must always be nonnegative, one can avoid the use of probability 
amplitudes in quantum mechanics. This proposal goes back to the work by Wigner who first introduced non- 

• positive pseudoprobabilities to represent quantum mechanics in phase space, and to the Moyal approach to quantum 
l/^ ■ mechanics ||^. 

• From a conceptual point of view, the elimination of the waves from quantum theory is in line with the procedure 
inaugurated by Einstein with the elimination of the aether in the theory of electromagnetism. 

The phase-space formulation of quantum mechanics provides a means of analyzing quantum-mechanical 

systems while still employing a classical framework. Moreover, a quantum mechanics without wave functions has 
been discussed in 

\ Recentely, the problem of quantum state measurement, initially posed by Pauli Q, received a lot of attention [p^ . 
^ ' The tomographic approach [pd] , p^ to the quantum state of a system has allowed to establish a map between the 
^ density operator (or any its representation) and a set of probability distributions, often called "marginals" . The latter 

have all the characteristics of classical probabilities; they are nonnegative, measurable, and normalized. 
^ Based on this connection, a classical-like description of quantum dynamics by means of "symplectic tomography" 

•'H . has been formulated jl^, providing a bridge between classical and quantum worlds. That is, the evolution of a 
r^J ' quantum system with continuous observables (namely, quadrature components of a field mode) was described in 
terms of a classical-like equation for a marginal distribution. Different aspects of this classical-like description using 
■ ■ ' tomographic probabilities were recently analyzed 

On the other hand, discrete observables, like spin or angular momentum, are so important in quantum mechanics as 
the continuous ones are. Hence, the tomography scheme for discrete variables was introduced [T^ , and the marginal 
distribution for rotated spin variables has been constructed JTzt , deriving an evolution equation for this function. 

Here, we would extend the approach by considering a spin-^ particle moving in a potential, then constructing 
the marginal distributions for space coordinates and spin projections and finally deriving the evolution equation for 
such probabilities, which would be the analog of the Pauli equation. It would also be a generalization of approaches 
attempted by us in previous papers [ p^ . 

Essentially, our aim is to eliminate the hybrid procedure of describing the dynamical evolution of a system, which 
consists of a first stage, where the theory provides a deterministic evolution of the wave function, followed by a hand- 
made construction of the physically meaningful probability distributions. If the probabilistic nature of the microscopic 
phenomena is fundamental, and not simply due to our ignorance, as in classical statistical mechanics, why should it 
be impossible to describe them in probabilistic terms from the very beginning? On the other hand, the language of 
probability, suitably adapted to take into account all the relevant constraints, seems to be the only language capable 
of expressing the fundamental role of "chance" in nature . 
The paper is organized as follow: 
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In Section 2, we review the general approach to construct known tomography schemes using density matrix in the 
specifically transformed reference frames. In Section 3, we derive the general evolution equation for tomographic 
probabilities (marginal distributions) which describe the quantum state instead of density matrix. In Section 4, the 
general scheme of tomography construction is used to re-derive the particular example of symplectic tomography, 
which is applied for measuring states depending on continuous quadrature. In Section 5, the general scheme is used 
to re-derive the construction of spin-state tomography. In Section 6, the general scheme of Section 2 is then applied 
to obtain tomographic probabilities in the combined situation described by spatial (multidimensional too) and spin 
variables. In Section 7, some examples are studied in the context of the probability representation of quantum 
mechanics. Section 8 concludes. 

We are using the natural unit {h = c = 1). 



II. GENERAL APPROACH TO QUANTUM TOMOGRAPHY 

In this section, we give a short review of the general principles used to construct a tomography scheme for measuring 
quantum states. Recently, we established JTof a quite general principle of constructing measurable probabilities, 
which determine completely the quantum state in the tomographic approach; more refined treatments then followed 
1^,^. Here, we apply our general approach to derive the evolution equation for the tomographic probabilities that 
is alternative in some sense to the Schrodinger equation for the wave function (or the quantum Liouville equation for 
the density matrix). 

Let us consider a quantum state described by the density operator p, which is a nonnegative hermitian operator, 
i.e., 

pt^p, Trp=l, (1) 

and 

{v\p\v)= p,,„ > 0. (2) 

We label the vector basis | v) in the space of pure quantum states by the multidimensional index v = {vi,V2, ■ ■ ■ ,vn) , 
where the number N shows the number of degrees of freedom of the system under consideration. Among indexes Vk , 
k = 1, . . . , N , there are continuous ones like position (or momentum) and discrete ones like spin projections. In this 
sense, the wave function ip{v) = (w | "0) of a pure state | ip) depends both on continuous and discrete observables. 
Formula (|^) can be rewritten by using the hermitian projection operator 

nv=\v){vi (3) 

in the following form 

p,,„ = Tr {n,p} . (4) 

The physical meaning of the projector IIi, is that it extracts the state | v) with given v (for example, with given 
position and spin projection), which is an eigenstate of the-commuting- hermitian operators V = {Vi, V2, • • ■ , Vn^ 

Vk\v)=vk\v). (5) 

In the space of states, there is a family of unitary transformation operators U{a) depending on the parameters a — 
(cti, . . . , (Tfc . . .), that can be sometimes identified with a group-representation operators. In these cases, the parameters 
a describe the group element. It was shown that known tomography schemes can be considered from the 

viewpoint of the group theory by using appropriate groups. More recentely this concept has been developed obtaining 
an elegant group theoretical approach to quantum state measurement Here, we formulate the tomographic 

approach in the following way. Let us introduce a "transformed density operator" 

p,^U-\a)pU{a). (6) 

Its diagonal elements are still nonnegative probabilities 

{Z \ Pa \ Z) = {{Z \ p \ Z)) EE W{Z, o) . (7) 
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Here, | z) is one of the possible vectors | w), while the symbol | z)) denotes the transformed vectors 



|z)> = C/(a)|z), (8) 
which in turn are eigenstates of thc-transformed-operators 

Z = l](a)Yij-^{o). (9) 
As a consequence of the unit trace of the density operator we also have the normalization condition 

dzw{z,a) = \. (10) 



Of course, in case of discrete indices, the integral in Eq. ( |10[ ) is replaced by a sum over discrete variables. 

Formula (^ can be interpreted as the probability density for the measurement of the observable Y in an ensemble 
of transformed reference frames labeled by the index cr, if the state p is given. Along with this interpretation, one can 
also consider the transformed projector 

n,(a)^C/(a)n,f7-i(cT)=|z))((z|, (11) 

the explicit expression for the probability w (z, cr) takes the form 

u;(z,a) =Tr{/5n4r7)} = Tr{/5 | |} . (12) 

These probability densities are also called "marginal" distributions as generalization of the concept introduced by 
Wigner . The tomography schemes are based on the possibility to find the inverse of Eq. ( |l^ ) . If it is possible to 
solve Eq. ([l^), considering the probability w {y,cr) as known function and the density matrix as unknown operator, 
the quantum state can be described by the positive probability instead of the density matrix. This property is the 
essense of state reconstruction techniques. In such cases, the inverse of Eq. ( p^ ) takes the form 

/5 = y w{z,a)K{z,a) dz da. (13) 

Thus, there exist a family of operators K (z, a) depending on both the variables z and parameters a such that the 
density operator is reconstructed, if the probability w (z, a) is known. It is worth remarking that transformations U {a) 
can form other algebraic constructions, which have no structure of groups [ p2[ . The only condition for the existence 
of a tomography scheme is the possibility to invert Eq. (p^). In the cases of optical tomography ||l^], symplectic 
tromography ||l3|, and spin tomography [ p7|j2^ , the sets of transformations U{a) have the structure of corresponding 
Lie groups (i.e., rotation, symplectic and spin). 



III. THE TIME EVOLUTION EQUATION 



We are now interested in finding the evolution equation for the probability w{z,a,t), in which t is the time 
parameter. Using Eq. (^2|) one has 

(14) 



^twiz,a,t) = TT^^[^tm]'nA^)} ■ 
On the other hand, the density operator satisfies the Liouville-Von Neumann equation 



(15) 

with H the system Hamiltonian. By inserting Eq.(|l5|) in (p^, and with the aid of Eq. (|l^) , we find the evolution 
equation for the probability w in a closed form 



dtw{z,a,t)= f dz' da' w{z',a',t)TT^^i K{z',a'),H n^(cr)| 



(16) 

Equation (|l^) represents the classical-like version of the Liouville-Von Neumann equation, thus, it would be the analog 
of the Pauli equation for a system with space and spin degrees of freedom. 
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IV. QUANTUM TOMOGRAPHY WITH CONTINUOUS VARIABLES 



Let us consider, in a one-dimensional system, an operator X as the linear combination of position q and momentum 
P 0,11 



X = + i/p, 



(17) 



which depends upon real parameters /i, v and, due to its hermiticity, is a measurable observable. Since the linear 
canonical transformation (17) belongs to the symplectic group Sp{2, R), the tomography scheme under discussion was 
called "symplectic tomography" j2^. 

The probability (marginal) related to the observable (^ is given by 

w{x,fi,i^) ^ {{x\ p\x)) , (18) 

where p is the system's density operator, while the eigenstates \ x)) of the operator ( p^ ) can be written as 

\x))^ J dq{q\x)) \q), (19) 

with I q) the position eigenkets. The wave function {q \ x)) can be easily calculated by using the following equality 

(q\X\x)) = {q\tiq + iyp\x)), (20) 
and then transforming it in a partial differential equation 

d 

x{q\x))^ixq{q\x))-iv—{q\x)). (21) 

The solution is 

{q\x)) = (^\'\^Ai^-q-l^qA . (22) 

It is worth noting that as soon as /x ^ 1 and — > 0, then | x)) — >| x) and the wavefunction ( p2| ) tends to S{q — x). 
Furthermore, Eq.(p!^) can be formally rewritten as 



w{x, fj., ly) =Tr ^ptl^{ti,iy)^ , 

where the transformed projector is given by 

U^ifi,!^) = U{p,iy)U^U-\p,iy) , fl, =| a;)(x | . 
Here, the transformation U (a) is chosen to be the symplectic group representation |lJ 

'iX 



U {p, v) = exp 



p^ q^ ^ 



exp 



[qp + pq) 



The rotation and scaling parameters (j) and A are related to p and v by the following formulas 

p = Xcos(j), j/ = A^^sin0, 



(23) 



(24) 



(25) 



(26) 



1 ■ , , 1 / 1 ± - Vi'^ 

- arcsm {2piy) , X = ±-]l 



This means that the marginal distribution w (x, /i, v) for this particular case of symplectic tomography is given by 
the relation 



w {x, p,v) = Tr < I x) (x I exp 



p q 



exp 



y {qp + pq) 













X pexp 






4)] 


exp 



(qp + pq) 



(27) 
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Such measurable probability can be explicitly expressed as 

w{x,^,i') = / dydk exp 



ikx H h iky^ 



piy + vk,y) , 



(28) 



where p{y + vk,y) = {y + i>k\p\y) is the representation of the density matrix over the position eigenkets. The marginal 
satisfy the following homogeneous property 



w{x, fi, Kv) — — w{x/ n, /i/k, J/) , 

w{x, KH, v) — — w{x/k, fj,, vj K) . 
K 

The above relation ( p^ ) can be inverted |2^] as 

P~ J dx dfjidv w{x^ iJi,v) K{x, jjL^v) , 
where the kernel operator takes the form 



(29) 



k{x,p,,u) = ^e^exp 



-ieX 



(30) 



(31) 



Here, e can be set equal 1; this freedom reflects the overcompleteness of information obtainable by means of all possible 
marginals (|2|) |2|,||. 

The multi-mode generalization j25j is straightforward, and the analog of formula (^^ holds with the following 
replacement 



I X) > I f ), X= {xi,X2, ■■■) 



p2 f 



^9 -2 

P i , Ti 



<t>2 



P2 , 12 



2 2 J " V2 2 J ^ \2 2 
A {qp + pq) > Ai {qipi + piqi) + A2 (92P2 + ^292) + 



(32) 



Relations of the parameters \k, (pk to the parameters iik, Vk are the same of Eq. (26) 



V. QUANTUM TOMOGRAPHY WITH DISCRETE VARIABLES 



Here, we consider a spin-j system. Following ||17|,^ we will derive the expression for the density matrix of a spin 
state in terms of measurable probability distributions. 

For arbitrary values of spin, let the spin state have the density matrix 



m = -j, + 1,. . . , j - 1, j, 



where 



and 



js I jm) = m I jm) , f \ jm) = j{j + 1) | jm) 



J J 

m=-j m' = -j 



(33) 
(34) 

(35) 



The operator p'^' is the density operator of the state under consideration. 

The general group construction of tomographic schemes was also used for spin tomography • The idea 

is to consider the diagonal elements of the density matrix p in another reference frame, i.e. rotated one. To this end 
we introduce a rotated measureble spin projection 



Ma, 13, 7) = Dia, l3, 7)j3^"' (a, /3, 7) , 



(36) 
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where the unitary rotation operator D depends on the Euler angles a, /3, 7 . The role of the observable Z is now played 
by the spin projection J3, while the rotation-transformation parameters are the Euler angles cti = a, <T2 = /3, = 7. 
The transformation U{a) is given by the matrix representation of the rotation group, i.e., the Wigner I^-function [ p7| . 
The marginals are 

w{s,a,f3,j) ^ {{js \ p \ js)) , (37) 

where the rotated spin states becomes 

\js))= ^ i?W*(a,/?,7) (38) 

Here the matrix elements -D,|;^/ ,„ {a, (3, 7) (Wigner /^-functions) are the matrix elements of the rotation-group repre- 
sentation 

7) = e*™'-^ e™" , (39) 

where 



41. = 



(j +m')!(j - m')!]^^^ ( . /^V" / . PY^ " p(m'-m,m'-(-m) 



(i + »™)!(j - "^)! 



cos^l sin^l F,^rj'"''"""^(cos/3), (40) 



with P^'^\x) the Jacobi polynomials [ pTt - 

Moreover, the transformed spin projector will be 

n«(a,/?,7) =^(a,/?,7) lisXjs | b-\a,Pn) ^]s)){{js\, (41) 

then, we have 

«;(s,a,A7)= ^ E ^KL(«,/3,7)pHU.i^E:(«,A7)- (42) 

Since 



^.:;;C„(«,/?,7) = (-1)" -™i^^^^.._™(a,/?,7) , (43) 
the marginal distribution really depends only on two angles, a and /?. Hence 

w{s,a,f3,j) ^ w{s,a,f3) , (44) 

which satisfies the normalization condition 

j 

^ «;(s,a,/3) = l. (45) 



As an example, for a spin-i state with spin projection +1/2, we have 



1 




(46) 



and the marginal distributions will be 

w (^s — i, a, P) = cos^ ^ , w (s = — i, a, /3) = sin^ ^ . (47) 



In Refs. |lj,^, in view of the properties of the Wigner ZJ-function and the Clebsch-Gordan coefficients, Eq. (42) was 
inverted and the density matrix was expressed in terms of the marginal distribution 
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2j js 



8^ 



where mi, m2 = —j, —j + 1, . . . , j and W^_^^^^^^^ are the Wigner-3 j symbols |27|| . The integration is performed over 
the rotation parameters, i.e. 

dn = da sin/?d/3 / dj . (49) 
Jo Jo JO 

Equation (^) can be presented in an invariant operator form . We systematically introduce the following notation, 
first for the function on the unit sphere 

ja 

™3 = -i3 

and then for the operator on the unit sphere 

if (a, P) = (2j3 + 1)' E El ^'"^1) '^J^^ ("' I . (51) 

mi = -j m2 = -j 

In order to write a final expression for the density operator, we introduce an operator on the unit sphere which 
contains a dependence on the measurable projection of the spin 

i:(^)(s,a,/3) = (-1)^ E W'i^fo^?'^ («,/?) • (52) 
Finally, we obtain a compact expression for the density operator, 

P'-'^ j ^ ^(^' ^^'^ (^' P)- (53) 

Formula ( |5^ ) admits of the following interpretation. To determine the spin state for a spin j , one has to experimentally 
measure the projection s of the spin for each direction specified by the angles a and (3, obtaining a distribution function 
w {s,a, (3). The sum on the r.h.s. of Eq.(p3|) for a given point on the unit sphere represents the average operator 
(i^T^j) (s,a,/3)). Then, the integral over the whole solid angle gives the desired density operator. Finally, we recognize 
that, for the spin case, the operator ( |52| ) plays the role of the operator K {z, a) of Eq. (|l3|), employed in the general 
scheme of Section 2. 



VI. THE GENERAL CASE 

We are now able to consider the case of a particle with — 1 spatial degrees of freedom, plus one spin-^ degree. 
In this case, the state vector | v) has the form 

I (f,™) =1 gi, . . .gjv_i)«) I i, m) , (54) 

where qis the eigenvalue of the position operator q and the spin projection m — (—1/2, 1/2) is the eigenvalue of the 
Pauli matrix a^. 

The transformation operator U{a) used to construct the tomography scheme, for this case, depends on 2{N — 1) 
parameters determining the symplectic transform, and on three Euler angles determining the spin rotation. 
The transformation operator U{a) of Eq. (^ becomes the product of operators 

Uia) = 0^-/^7 (Mfc, i^k) ® U (a, (3, 7) . (55) 
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For the case of spin-i, the representation of the rotation group is given by 



7) 

which determines the operator 

[7(a,/3,7)- 



e*"/2 cos {(5/2) e-^'^l'^ -e^'^/^ (-^/2) e*^/^ 
e*"/2 sin (/3/2) e^"/^ cos (/5/2) e^^'"^ 



1/2 1/2 

E E ^^/l(«,/3,7)li"^i)(im2 

mi=-l/2 m2 = -l/2 



(56) 



(57) 



The marginal distribution w(z^a) ( p^ ) will depend on — 1 continuous (noncompact) variables z\ — xi, 
zn-i = xat-i, and one discrete spin projection z^r = s, as well as on parameters /ife, and on Euler angles a, /?. 
The dependence of the marginal distribution on the Euler angle 7 disappears, as it was shown in the previous section, 
due to the structure of Wigner 13- functions. 

In order to get an analog of the Pauli evolution equation for the marginal distribution, we consider the general 
equation (p^ where the operator K {z' , a') has the form 



k{z\a') = 



(58) 



Here, the operator K{xk, fJ,k, i^k) has the form of Eq.(^l|) with e — 1, and the operator (s,a, (3) is given by 

formula (|5^ ) with j = 1/2. Moreover, we have to introduce the marginal distribution w {x, jl, v, s, a, /3, t) describing a 
state of spin-i particle which depends on the continuous variables x, discrete spin projection s, symplectic reference 
frame's labels pi and v, and Euler angles a and [3. Then, for a given Hamiltonian H the general equation (^6|) takes 
the form of a Pauli-like equation equation 



dtw{x, fl,iy,s,a,f3,t) 



dX'dfl'dv'dVL' w{x', p', v', s', a', /?', t) 



1/2 

- E 

s' = -l/2 ' 

X Q{x, fi, V, s, a, /3; x' , 0, v' , s', a' , (3') , 



where 



e 



((^,^1 \®^~,'k{x'^.n',y'u)®k'^''^H-^'.c.\ fi'),H\ \x,s)) 



(59) 



(60) 



The structure of the derived Pauli-like equation for probability distributions depends on the particular tomography 
schemes we have considered. Obviously, it would be useful to find the schemes which give the simplest form for such 
dynamical equ ation, nevertheless this is a nontrivial problem related to the possibility of finding properly transformed 
projector (ll|). The latter are investigated in Ref. |20|, but for different purposes. 



A. Limit cases 

We want now to consider two limiting cases of the above equation (|5£ 
First of all we consider the (one-dimensional) spatial case only with free motion 

i? = y- (61) 
The spin part does not contribute since H does not contain the spin operators, that is 

«;(.', a', /3',-)((s I («',«',/?') \ s)) ^ ^ ^ (a, /3, 7)i^?> K A 7) 

mi,m2 = -j 

E E E (-)"""^'(2j3 + i)^^v^/i,ioW^//:^™.™3 

X / _^„(s',a',/3',-)i?(tU«''/^''V)=«^(^^,«,/3,-), (62) 
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where — indicates other possible variables. Then, for what concerns the spatial part, it is important to calculate the 
commutator between the kernel and the Hamiltonian, given by 



iu' a if' p •:-2 



(63) 



Now, one can write 



dtwix,iJ,,i',t) = ^ / dx'dii'diy'wiX',ii',i'',t)e-'^'+"'''''/^ 

4:71 J 

X J dq{{x I e*^'«e-'^^|g)(g|(-2M'p - m'') | x)) (64) 

By using the explicit form for the wave functions (q \ x)) (p2|), toghether with the homogeneous property (p9|), it is 
possible to reduce the above equation to a very simple form 



dtw — ^d^w 

which was derived in a different way in Ref. jT^ . 

As a second case we study the dynamics of spin-^ degree only. The Hamiltonian we wish to consider is 



H = 



a 
c 



(65) 



(66) 



Of course, the spatial degree is not affected, so its variables can be disregarded; this also results from the fact that 

-!- / dx'd^i'di^'wix',fi',iy',^){{x\e"'''^e"''P \ x))e~"'+'^''''' /'^ = w{x, (67) 
27r J 

In this case the relation between the transformed spin state projection and the untransformed one is given by 

Again, the central task is the calculation of the commutator between the kernel and the Hamiltonian. It is easy to 
see that 



{{s\ k<^'I^Hs',a\nH \s)) = {-l)-'' Y^Wl; 



1/2 1/2 J3 



-s' 



(2^3 



33=0 



1/2 



J 3 



X E ^-r^ E ^^tl 7') <C 

mi5^m2,-l/2 ™3 = -j3 

X „ e)i?( ^)^(l/2) * ^) 



1/2 1/2 i3 
m2 m3 



(69) 



Due to the properties of the Wigner-3j symbols we may see that the terms with j3 = 0, 1, and = do not give 
contributions; moreover, changing the value of s', it changes only the sign. Thus, we will get 



dt'w{\,a,l3,t) 



dn' 



X -(a — c) sin/3' sin/3sin(a — a') 



and by using the normalization condition it can be rewritten as 

dfl' 
8^ 



dt 'w{s,a, (3,t) ~ 3{a — c) sin P j ^ w(s, a', /?', i) sin/3' sin(a — a') 



(70) 



(71) 



which is similar to that derived in Ref. |17| (the differencies are due to the degeneracy of the spin-i systems). It 
should be noted in the above equation that the argument s is the same in both sides; this is consistent with the fact 
that H in Eg. (|66|) does not mix states with different s. On the other hand it can be easily checked that the sum over 
s at r.h.s. of Eq.(|7l]) causes the integral to become zero; this is consistent with the fact that at l.h.s. we will obtain 
the time derivative of a constant. Also, if a = c, the r.h.s. of Eq.(^l]) will be zero since the Hamiltonian ( |66| ) will be 
proportional to the identity and will not produce any evolution. 
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VII. EXAMPLES 



In the previous section, we discussed the probabihty for the joint measurement of the spin and spatial variables. 
Therefore, here we would like to consider some examples involving both variables. 
At first we consider a system with the following hamiltonian 



H-\{f+e) + m\\-\-\){-\\) 



(72) 



It could describe e.g. one vibrational degree of a trapped electron plus its spin |28|. The measurability of marginals in 
this system is investigated in Ref. |E9|. Here, as a straigthforward extension of the arguments of Sec. VI. 1 we obtain 



dt w{x, /i, V, s, a, P) — {^idi, — i^df_i) w{x, /i, v, s, a, (3) 

8^ 



(73) 



+ 6sin/3 / ^ wlx, fj,, v, s, a' , (3' ,t) sin P' sm{a — a') . 



Let us now consider an initial entangled state like 

1 



*(0) = ^(|0)®|-i)+|l)®|i)). 



(74) 



where | n) represents the number eigenstate of a harmonic oscillator. At Eq.([74[) corresponds the following marginal 



where 



w{x, V, s,a,P,t ^ 0) = - [woaii + twiin + wout + wionl > 



1 



wooii 

Willi 
WOUl 



exp 



2x^ 

— , cxp 

^y^^ip'^ + 1/2)3 

^/2x(i> — in) 
i — — — exp 



/i2 + jy2 



^2 _|_ 



^7r(^2 ^ ^2)3 

Then, the solution of the Pauli equation (|73|) is 



Ij2 ^ jy2 



i?(^5^*(a,/3,7)i^^'5^(a,/3,7), 

2 2 

"2 ^2 

* 2 ^2 



(75) 

(76) 

(77) 

(78) 
(79) 



w{x, ^, s, a, p,t) ^ - [wooii + ti'mi + woiiie^'* + wwu 



-3itl 



(80) 



As a second example, we want to consider the case of Landau levels ||30|| , i.e. a charged particle moving in a classical 
magnetic field B being time-independent and axial symmetric. The particle's movement along the axis being free, 
instead the Hamiltonian of the transverse motion reads 



H = 



1 



(pl - il) + (p2 - 



A = 



B X 



(81) 



where f— (gi, 172) is the radius- vector of the particle's center, pi and p2 are the particle's momentum components in 
the transverse plane. Having B along the third axis and choosing | _B |= 2, we get 



H 



{Pl+Pl + 9? + €) + {Pl<l2-P2qi) + - \-\){-\\) 



(82) 



In this case the kernel Q of Eq. (pOh is given by 
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47r2 J J V2 

X \ - [^iXijvi + v[qi - [l!^ jl - v'^ 12) - {X2 - Ai292)K/'^2 

1 1=1 

+ (M29i - ^21^1) + (2^1 - Mi9i)^'2/^^i - (m'i92 - \ Ss.s' S{Q - Q') 



6 
8^ 

As nontrivial example we also consider here an initial state which is the entangled superposition 



+ ^— jsin/3sin/3'sin(a — a')6{x — x')5{fl — pi')5{v — if) . (83) 



f(0) = i=[|00)®Ki)+|10)®|i)] . (84) 

It leads to nonfactorisablc marginal 

w (f , /I, V, s, a, /3, t = 0) = i [woooou + WioioTT + ^"ooioiT + fi'ioooTil > (85) 

where 

Wm n2 n\ n'^ ,mi ma = "'ni n'^ X D^^^/^^ *(a, 7)^i Pi l) ■ (86) 

Here, mi = — ^ (mi = i) replaces downarrow (uparrow) while the spatial part w„j^„2n'j^n^ is explicitely calculated in 
the Appendix. 

It is now easy to see that the solution of the Eq.(^9|) with the kernel (|8^), subject to the above initial condition, is 
w (f , fl, u, s, a, (3,t)^^ [wqooou + wioiott + wonini^e^'* + wiooone"^**] ■ (87) 



VIII. CONCLUSION 

We conclude that it is possible to obtain an evolution equation for the tomographic probabilities (marginal distri- 
butions) of an arbitrary tomography scheme. The main result of our paper is the analog of the Pauli equation for 
spin-i particle. 

The explicit expression for the marginal distribution for a trapped particle as well as for Landau levels has been 
studied. It results that in the nonstationary case they obey the analog of the Pauli equation. 

The examples considered demonstrate that the usual problems of conventional quantum mechanics can be cast into 
the form in which only positive probabilities are used to describe quantum states and their evolution. A possible 
disadvantage of the approach proposed is a complicated evolution equation (|5|), but, perhaps, this is the price one 
ought to pay for the possibility of describing quantum objects in terms of classical probabilities. 

Anyway, our argumentations can constitute a step further from the Bohr position |]3l|| about the inapplicability 
of classical modes of description in the quantum domain. In fact, while we belive that quantum mechanics is not 
classical physics in disguise, we retain (some) classical concepts still applicable against counterintuitive notions like 
complex statefunctions. 

We also belive that the developed classical-like formalism could be applied to describe quantum mechanical para- 
doxes, because usually, if there is a paradox in quantum mechanics, there should also be a classical one, perhaps worse 
p2[ . These aspects will be investigated in a forthcoming paper as well as the extension of the presented approach to 
the relativistic domain [3^ ], in order to find an analog of the Dirac equation. 
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APPENDIX 



The wave function of the particle's coherent state in a magnetic field B is [Q] 

^a,p{qi,q2) = ■^'^^P I ^ 2 ^ 2 2 + + *92) + ia (qi - iq2) 

where qi and q2 are the particle's coordinates and a and f3 complex numbers. 
The coherent state ( p8| ) is the superposition of number states [|4j 

(.1, = exp (^^^ E — — • m 

In view of the general relationship between the marginal distribution and wave function , we have 

1 



47r2 I viV2 



, -iyffj-i iyixi , iyj^j-2 iy2X2\ , \ ^ ^ 
cxp — h — *a/3 (yi,y2) dyi dy2 

2vi Vi 2V2 V2 ) 



(90) 



where parameters /ii, /i2, i^2, as usually, mark reference frames; then, one obtains for the marginal distribution of 
the particle's coherent state without spin in a magnetic field the following expression 

exp[-|a|2-|/3|2_z(«/3_a*;3*)] 

[Xi,X2,\i.\,VY,[l2,V2) = — 



X exp 



y 2VI + 1/2) 

^ (i/2 + 1^2) (a«^2 + j/3i^2 - •1x2)'^ + (;^2 - 'ifJ-2) {a*V2 - if3*i^2 + 1x2)^ 1 
Multiplying (|9l| ) by exp (jap + and expanding the expression obtained into the power series, we arrive at 



(91) 



WafS {xi,X2, 1^1,^2, 1^2) e^"^ e'"! = > ======= —. (92) 

■~ , „ ^/nT\n2\n\lnJ. 

ni ,712 ,n2— ^ i z 

Taking into account the property of the generating function for multivariate Hermite polynomials js^, namely, 

exp|-i.^M.^+.^AfcU E ^ /-a 

ni,n2,n'^ ,712—0 

where the vector u has components u — {a, a* , /3, /?*), and comparing (^) with (^3|), we obtain 

1 



(^711 772 7i; 77;^ {X1,X2, fJ-l,!^!, fJ-2, V2) = 



[Ml /^A 

11 712 n'^ ' '"■2 \ / 

+ ^2 + ^^2 7 v/"i!"2!ni!n2! 



^2 X 



where the 4x4 matrix M reads 



The 2x2 matrices M''') are given by 



12 



1 Vj+ii-Yt^i 



I. ) "fc+i,even ; ' — ^t'J i 



(96) 



i=i 



(_)'+0-l)('-/2-l) ^_ (-)'-i 



<5fc+i,ovcn , J' = 2, 4, 



The argument of the multivariate Hermite polynomials C = (Cij Ci j C27 C2) expressed in terms of the parameters as 
follows 



Ci = / o cxp I tan — 



C2 



iX2 



2:2 



^2/ VmT 



exp I tan — 



cxp I i tan ^ — ) , 



exp I tan — 

1^2 



(97) 



Taking ni — 712 and n'l ~ 712 wc obtain the marginal distribution i(;„„'(a;i,a;2,/ii,J^i,/i2,f2) for the Landau level states 
I nn') 



Wnn' iXl,X2,^J,l,Vl,^^2,V2) = W„„„'„' (xi , ^2 , Ml 7 ^^1 , ^2 , ^^2) 



exp 



H, 



{M} 



..(c) 



^ VR+M?)(i?+Ml)n! n'! 
where n is the main quantum number and n' ~ n = I is the angular momentum quantum number. 
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